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ABSTRACT 

Recent advances in understanding of magnetohydrodynamic (MHD) turbu- 
lence call for substantial revisions in the picture of cosmic ray transport. In this 
paper we use recently obtained scaling laws for MHD modes to calculate the 
scattering frequency for cosmic rays in the ISM. We consider gyroresonance with 
MHD modes (Alfvenic, slow and fast) and transit-time damping (TTD) by fast 
modes. We provide calculations of cosmic ray scattering for various phases of in- 
terstellar medium with realistic interstellar turbulence driving that is consistent 
with the velocity dispersions observed in diffuse gas. We account for the turbu- 
lence cutoff arising from both collisional and collisionless damping. We obtain 
analytical expressions for diffusion coefficients that enter Fokker-Planck equation 
describing cosmic ray evolution. We calculate the scattering rate and parallel 
spatial diffusion coefficients of cosmic rays for both Alfvenic and fast modes. We 
conclude that fast modes provides the dominant contribution to cosmic ray scat- 
tering for the typical interstellar conditions in spite of the fact that fast modes 
are subjected to damping. We show that the efficiency of the scattering depends 
on the plasma j3 since it determines the damping of the fast modes. We also 
show that the streaming instability is modified in the presence of turbulence. 

Subject headings: acceleration of particles-cosmic rays-ISM: magnetic fields- 
MHD-scattering-turbulence 

1. Introduction 

Most astrophysical systems, e.g. accretion disks, stellar winds, the interstellar medium 
(ISM) and intercluster medium are turbulent with an embedded magnetic field that influ- 
ences almost all of their properties. High conductivity of the astrophysical fluids makes the 
magnetic flelds "frozen in", and influence fluid motions. The coupled motion of magnetic 
fleld and conducting fluid holds the key to many astrophysical processes. 
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The propagation of cosmic rays (CRs) is affected by tfieir interaction witfi magnetic 
field. This field is turbulent and therefore, the resonant interaction of cosmic rays with 
MHD turbulence has been discussed by many authors as the principal mechanism to scatter 
and isotropize cosmic rays (Schlickeiser 2002). Although cosmic ray diffusion can happen 
while cosmic rays follow wandering magnetic fields (Jokipii 1966), the acceleration of cosmic 
rays requires efficient scattering. For instance, scattering of cosmic rays back into the shock 
is a vital component of the first order Fermi acceleration (see Longair 1997). 

While most investigations arc restricted to Alfven modes propagating along an external 
magnetic field (the so-called slab model of Alfvcnic turbulence), obliquely propagating MHD 
modes have been included in Fisk et al. (1974) and later studies (Bieber et al. 1988, Pryadko 
& Petrosian 2000). The problem, however, is that the Alfvenic turbulence considered in their 
studies is isotropic turbulence, and this contradicts to the modern understanding of MHD 
turbulence (Goldreich & Shridhar 1995, see Cho, Lazarian & Yan 2002 for a review and 
references therein). 

A recent study (Lerche & Schhckeiser 2001) found a strong dependence of scattering 
on turbulence anisotropy. Therefore the calculations of CR scattering must be done using 
a realistic MHD turbulence model. An important attempt in this direction was carried out 
in Chandran (2000). However, only incompressible motions were considered. On the con- 
trary, ISM is highly compressible. Compressible MHD turbulence has been studied recently 
(see review by Cho & Lazarian 2003a and references therein). Schlickeiser & Miller (1998) 
addressed the scattering by fast modes. But they did not consider the damping, which is es- 
sential for fast modes. In this paper we discuss in detail the various damping processes which 
can affect the fast modes. To characterize the turbulence we use the statistics of Alfvenic 
modes obtained in Cho, Lazarian & Vishniac (2002, henceforth CLV02) and compressible 
modes obtained in Cho & Lazarian (2002, henceforth CL02, 2003b,c). 

Yan & Lazarian (2002) used recent advances in understanding of MHD turbulence 
(CL02) to describe cosmic ray propagation in the galactic halo. In this paper we under- 
take a comprehensive study of cosmic ray scattering rates. In §2, we describe the statistics 
of Alfvenic and compressible turbulence in various conditions, including both high (5 and 
low j3 cases. In §3, we describe the resonant interactions between the MHD modes and CRs. 
The scattering by Alfven modes and fast modes is presented in §4. We apply our results to 
different phases of ISM, including Galactic halo, hot ionized medium (HIM), warm ionized 
medium (WIM) and also partially ionized medium. In §5, we describe the streaming insta- 
bility in the presence of background turbulence. Discussion of our results is provided in §6 
while the summary is given in §7. 



-3- 



2. MHD cascade and its damping 

2.1. MHD turbulence cascade 

A substantial shortcoming of earlier studies was that the crucial element for cosmic ray 
scattering, namely, MHD turbulence model was taken somewhat ad hoc. The goal of the 
current paper is to use the recent advances in quantitative description of MHD turbulence 
to quantify scattering of cosmic rays. 

MHD perturbations can be decomposed into Alfvenic, slow and fast waves with well- 
known dispersion relations (see Alfven & Falthmmar 1963). Alfvenic turbulence is considered 
by many authors as the default model of interstellar magnetic turbulence. This is partially 
motivated by the fact that unlike compressible modes, the Alfven ones are essentially free 
of damping in fully ionized medium (see Ginzburg 1961, Kulsrud & Pearce 1969). Impor- 
tant questions arise. Can the MHD perturbations that characterize turbulence be separated 
into distinct modes? Can the linear modes be used for this purpose? The separation into 
Alfven and pseudo- Alfven modes is the cornerstone of the Goldreich-Sridhar (1995, hence- 
forth GS95) model of turbulence. This model and the legitimacy of the separation were 
tested successfully with numerical simulations (Gho & Vishniac 2000; Maron & Goldreich 
2001, GLV02). Separation of MHD perturbations in compressible media into fast, slow and 
Alfven modes is discussed in GS95, Lithwick & Goldreich 2001, GL02). The actual decom- 
position of MHD turbulence into Alfven, slow and fast modes was performed in CL02, Cho 
& Lazarian (2003, henceforth CL03), who also quantified the intensity of the interaction 
between difi^erent modes (see below). 

Turbulence of Alfven modes is the default model of MHD turbulence for many resources. 
Other models are considered less important due to damping. However, within an earlier re- 
search it was frequently forgotten that, unlike hydrodynamic turbulence, Alfvenic turbulence 
is anisotropic, with eddies elongated along the magnetic field (see Higdon 1984, Shebalin et 
al 1983). This happens because it is easier to mix the magnetic field lines perpendicular to 
the direction of the magnetic field rather than to bend them. The GS95 model describes 
incompressible Alfvenic turbulence, which formally means that plasma (3 — Pgas/Pmag, the 
ratio of gas pressure to magnetic pressure is infinitely large. The corresponding scaling can 
be easily obtained. For instance, calculations in CLV02 prove that motions perpendicular 
to magnetic field lines are essentially hydrodynamic. As the result, energy transfer rate 
due to those motions is constant Ek ~ ^l/rk, where is the energy eddy turnover time 
~ {vkk±)~^, where k± is the perpendicular component of the wave vector k. The mixing 
motions couple to the wave-like motions parallel to magnetic field giving a critical balance 
condition, i.e., k±Vk ~ k\\VA, where k\\ is the parallel component of the wave vector k, Va 



-4- 



is the Alfven speed^. Prom these arguments, the scale dependent anisotropy oc k_[ and 
a Kolmogorov-Ukc spectrum for the perpendicular motions V}^ oc k~^/^ can be obtained (see 
Lazarian & Vishniac 1999). 

It was conjectured in Lithwick & Goldreich (2001) that GS95 scahng should be approxi- 
mately true for Alfven and slow modes in moderately compressible plasma. For magnetically 
dominated, the so-called low j3 plasma, CL02 showed that the coupling of Alfvenic and com- 
pressible modes is weak and that the Alfvenic and slow modes follow the GS95 spectrum. 
This is consistent with the analysis of HI velocity statistics (Lazarian & Pogosyan 2000, 
Stanimirovic & Lazarian 2001) as well as with the electron density statistics (see Arm- 
strong, Rickett & Spangler 1995). Calculations in CL03 demonstrated that fast modes are 
marginally affected by Alfven modes and follow acoustic cascade in both high and low (5 
medium. In what follows, we consider both Alfven modes and compressible modes and use 
the description of those modes obtained in CL02, CL03 to study CR scattering by MHD 
turbulence. 

The distribution of energy between compressible and incompressible modes depends, in 
general, on the way turbulence is driven. CL02 and CL03 studied generation of compressible 
perturbations using random incompressible driving. As the result we obtained an expression 
that relates the energy in fast ~ 5Vf and Alfven ~ SV^ modes, 

{5Vfl5VAf ~ 5Va X VaI{VI + C|), (1) 

where Cg is the sound speed. This relation testifies that at large scales incompressible 
driving can transfer an appreciable part of energy into fast modes. However, at smaller 
scales the drain of energy from Alfven to fast modes is marginal. Therefore the cascades 
evolve without much of cross talk. Naturally a more systematic study of different types of 
driving is required. In the absence of this, in what follows we assume that equal amounts of 
energy are transfered into fast and Alfven modes when driving is at large scales. 

At small scales turbulence spectrum is truncated by damping. Various processes can 
damp the MHD motions (see Appendix A for details). In partially ionized plasma, the 
ion-neutral collisions are the dominant damping process. In fully ionized plasma, there are 
basically two kinds of damping: collisional or collisionless damping. Their relative impor- 
tance depends on the mean free path in the medium (Braginskii 1965), 

Imfp = VthT = 6 X 10^^ {T/mmKf/n. (2) 



note that the linear dispersion relation is used for Alfven modes. 
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If the wavelength is larger than the mean free path, viscous damping dominates. If, on 
the other hand, the wavelength is smaller than mean free path, then the plasma is in the 
coUisionless regime and coUisionless damping is dominant. 

To obtain the truncation scale, the damping time F^^ should be compared to the cas- 
cading time Tfe. As we mentioned earlier, the Alfvenic turbulence cascades over one eddy 
turn over time {k±Vk)~^ ~ {k\\VA)~^- The cascade of fast modes takes a bit longer: 

r, = co/k'vl = (k/L)-'/' X (3) 

where V is the turbulence velocity at the injection scale, Vph is is the phase speed of fast 
modes and equal to Alfven and sound velocity for high and low /3 plasma, respectively 
(CL02). If the damping is faster than the cascade, the turbulence is truncated. Otherwise, 
for the sake of simplicity, we ignore the damping and assume that the turbulence cascade 
is unaffected. As the transfer of energy between Alfven, slow and fast modes of MHD 
turbulence is suppressed, we consider different components of MHD cascade independently. 

We get the cutoff scale kc by equating the damping rate and cascading rate TkT^ ~ 1. 
Then we check whether it is self-consistent by comparing the kc with the relevant scales, e.g., 
injection scale, mean free path and the proton gyro-scale. 

Damping is, in general, anisotropic, i.e., the linear damping (see Appendix A) depends 
on the angle between the wavevector k and local direction of magnetic field B. In the typical 
ISM this happens at the sufficiently small scales where the direction of B is well defined (see 
discussion in the review by Cho, Lazarian & Vishniac 2003). 

Note that unless randomization of wave vectors is comparable to the cascading rate 
the damping scale gets angle-dependent. Consider fast modes, that will be shown to be the 
most important for cosmic ray scattering. Their non-linear cascading can be characterized by 
interacting wavevectors k that are nearly collinear (see review by Cho, Lazarian & Vishniac 
2003). The possible transversal deviation 6k can be estimated from the uncertainty condition 
^'^tcas ~ 1; where 5uj ~ Vph5k{5k/k). As teas ~ l/fc^'ffc ^ k~^l'^ (see CL02) the bk scales 
as k^l'^. Therefore the deviation in k direction in the course of interactions decreases with 
the increase of the wavenumber as Skjk ~ k~^l'^^ and the number of interactions required 
to randomize k increases as k}!'^ with the increase of the wavenumber. This means that the 
randomization of k gets marginal for large k. In the presence of anisotropic damping this 
results in anisotropic distribution of fast mode energy at small k. 

With this input at hand, it is possible to determine the turbulence damping scales for a 
given medium. For modes in fully ionized medium one should compare the wavelength and 
mean free path before determining which damping we should apply. However, we are dealing 
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with the turbulence cascade. Therefore, the situation is more comphcated: If the mean 
free path is larger than the turbulence injection scale, we can simply apply the collisionless 
damping to the whole inertial range of turbulence. In general one should compare the viscous 
cut off and the mean free path. If the cutoff scale is larger than the mean free path, it shows 
that the turbulence is indeed cut off by the viscous damping. Otherwise, we neglect the 
viscous damping and just apply the collisionless damping to the turbulence below the mean 
free path^. By comparing different damping, we find the dominant damping processes for 
the idealized ISM phases (see tablel). 

2.2. Statistics of fluctuations 

Within random-phase approximation, the correlation tensor in Fourier space is (see 
Schlickeiser & Achatz 1993) 

< B,{k)B*{k') > /Bl = 5(k - k')M,,(k), 
< v,{\i)B*{k!) > /VaBo = 6{k - k')Q,(k), 

< ^.(k)^;(k') > /V2 = 6{k - k')K,,{k), (4) 

where B^^p, Va,i3 are respectively the magnetic and velocity perturbation associated with 
the turbulence. For the balanced cascade we consider, i.e., equal intensity of forward and 



^We shall show that for some angles between B and k the damping may result in turbulence transferring 
to collisionless regime only over a limit range of angles. 



ISM 


galactic halo 


HIM 


WIM 


WNM 


CNM 


MC 


T(K) 
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lmfp{cm) 


4 X 10^9 
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6 X 10^2 


8 X 10^1 


3 X 10^ 


10^ 


L(pc) 


100 


100 


50 


50 


50 


50 




5 


2 


5 


5 


5 


15 


damping 


collisionless 


collisionless 


coUisional 


neutral-ion 


neutral-ion 


neutral-ion 



Table 1: The parameters of idealized ISM phases and relevant damping. The dominant 
damping mechanism for turbulence is given in the last line. HIM=hot ionized medium, 
CNM=cold neutral medium, WNM=warm neutral medium, WIM=warm ionized medium, 
DC=dark cloud. 
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backward modes, Cij(k) = 0. 

The isotropic tensor usually used in the literature is 

Mij{k) = Co{5ij - kikj/e}k-''/\ (5) 

The normalization constant Cq here can be obtained if the energy input at the scale L 
is defined. Assuming equipartition Uk = J dk^ Yli=i ^u^q/^^ ~ Bq/Stt, we get Cq = 
L~^/^/127r. The normahzation for the following tensors are obtained in the same way 

The analytical fit to the anisotropic tensor for Alfven modes, obtained in CLV02 is. 



M,,(k) 



-1/3 



Gtt 



I,jkZ''^'eM-L'^%\\/kT), 



(6) 



where lij 



{Sij — kikj/k"^} is a 2D tensor in x — y plane which is perpendicular to the 
magnetic field, L is the injection scale, V is the velocity at the injection scale. Slow modes 
are passive and similar to Alfven modes. 

According to CL02, fast modes arc isotropic and have one dimensional energy spectrum 
E{k) (X A;~^/^. In low (3 medium, the corresponding correlation is (YL02) 



Kij{k) 



Stt ^ 



cos^ 6 
1 



(7) 



where 9 is the angle between k and B, Hij = kikj/k^ is also a 2D tensor va. x — y plane. The 
factor cos^ 6 represents the projection as magnetic perturbation is perpendicular to k. This 
tensor is different from that in Schlickeiser & Miller (1998)^. For isotropic turbulence, the 
tensor of the form oc Ek{5ij — kikj/k^) was obtained to satisfy the condition k • 5B = (see 
Schlickeiser 2002). Nevertheless, the fact that S'B in fast modes is in the k-B plane place 
another constraint on the tensor so that the term Sij doesn't exist. 

In high (3 medium, fast modes in this regime are essentially sound modes compressing 
magnetic field (GS95, Lithwick & Goldreich 2001, CL03). The compression of magnetic field 
depends on plasma /3. The corresponding x-y components of the tensors are 



M,,(k) 
KM 



7- -1/2 

27r ^ 



cos2 e/p^ 

1/P 



(8) 



^Here we give only the x,y component for the perturbation, the solenoidal condition wiU be satisfied if 
the z component is added. 
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The velocity perturbation in high /3 medium is radial, i.e., along k, thus we have the factor 
sin^ ^ and also C^/V^ — 2/j3 from the magnetic frozen condition a;5B ~ k x {y^ x B). In 
high [3 medium, the energy of fast modes are reduced by V^/C|, this gives another 2/(3. We 
use these statistics to calculate cosmic ray scattering arising from MHD turbulence. 



3. Interactions between turbulence and particles 

Particles get into resonance with MHD perturbations if the resonant condition is fulfilled, 
namely, uj — k\]Vfi + riQ {n = ±1,2...), where ui is the wave frequency, fl = flo/j is the 
gyrofrequency of relativistic particle, /i — cos^, where ^ is the pitch angle of particles. 

Basically there are two main types of resonant interactions: gyroresonance acceleration 
and transit acceleration (henceforth TTD). The latter requires longitudinal motions and it 
only operates with compressible modes. Gyroresonance occurs when the Doppler shifted 
frequency of the wave in the particle's guiding center rest frame Ugc = uo — k\\v^ is a multiple 
of the particle gyrofrequency, and the rotating direction of wave electric vector is the same 
with the direction for Larmor gyration of particle. For high energy particles, the resonance 
happens for both positive and negative n. TTD happens due to the resonant interaction 
with parallel magnetic mirror force —{mv\/2B)V\\Q. For small amplitude modes, particles 
should be in phase with the wave so as to have a secular interaction with wave (Schlickeiser 
& Miller 1998). This gives the Cherenkov resonant condition u — k\\v\\ — 0. 

We employ quasi-hnear theory (QLT) to obtain our estimates. QLT has been proved 
to be a useful tool in spite of its intrinsic hmitations (Schlickeiser & Miller 1998, Miller 
1997). For moderate energy cosmic rays, the corresponding resonant scales are much smaller 
than the injection scale. Therefore the fluctuation on the resonant scale 5B <C Bq even if 
they are comparable at the injection scale. QLT disregards diffusion of cosmic rays that 
follow wandering magnetic field lines (Jokipii 1966) and this diffusion should be accounted 
separately. If mean magnetic field is larger than the fiuctuations at the injection scale, we may 
say that the QLT treatment we employ defines parallel diffusion. Obtained by applying the 
QLT to the coUisionless Boltzmann-Vlasov equation, the Fokker-Planck equation is generally 
used to describe the evolvement of the gyrophase-averaged distribution function /, 



dt d^i V ""dn 



D ^V-- 
dp J dp 



p 



D ^ 



D ^ 



where p is the particle momentum. The Fokker-Planck coefficients D^^, D^p, Dpp are the fun- 
damental physical parameter for measuring the stochastic interactions, which are determined 
by the electromagnetic fluctuations (Schhckeiser 1993): 
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Adopting the approach in (Schhckeiser 1993) and taking into account only the dominant 
interaction at n = ±1, we can get the Fokker-Planck coefficients (see Appendix B), 



D 

n 



pp 





dk^?>{kf)\\ - a; + ±Q) 



2J.(^)Jo 



+ e 



-i2<f> 



(9) 



where ^ = 1 for kmin = L~^, kmax = ^o/i^th corresponds to the dissipation scale, m — ^mn 
is the relativistic mass of the proton, v±_ is the particle's velocity component perpendicular 
to Bq, = arctan(A;j^/A;a.), C^TZ — {x ^ iy)/\/2 represent left and right hand polarization. 

The delta function — uj + nVl) approximation to real interaction is true when 

magnetic perturbations can be considered static^ (Schhckeiser 1993). For cosmic rays, k\\v\\ ^ 
uj so that the resonant condition is just k\]V^-\- nil = 0. Prom this resonance condition, we 
know that the most important interaction occurs at k\\ — k\\^res — ^/v\\. This is generally 
true except for small ji (or scattering near 90°) . In this case, these particles move at a speed 
close to Alfven speed along the magnetic field so that we have to take into account the 
dynamics of the turbulence (Schhckeiser 1993, Bieber, Matthaeus & Smith 1994, Yan & 
Lazarian 2003). 



4. Scattering of cosmic rays 

4.1. Scattering by Alfvenic turbulence 

As we discussed in §2, Alfven modes arc anisotropic, eddies arc elongated along the 
magnetic field, i.e., k±_ > k\\. The scattering of CRs by Alfven modes is suppressed first 
because most turbulent energy goes to k± due to the anisotropy of the Alfvenic turbu- 
lence so that there is much less energy left in the resonance point k\\^res — ^/'^^|| ~ '"Z^- 



^Cosmic rays have such high velocities that the slow variation of the magnetic field with time can be 
neglected. 
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Furthermore, k± ^ A;|| means k± ^ r^^ so that cosmic ray particles have to be inter- 
acting with lots of eddies in one gyro period. This random walk substantially decreases 
the scattering efficiency. Mathematically, this effect is embodied by the Bessel functions in 
the Fokker-Plank coefficients. We know that the asymptotic expression of Bessel function 
Jn{x) ~ Xj \fx. So the suppression will be of the order (f2/A;_L?;^) 2 = {k±rL)'2. Noticing that 
k±v±/fl = /c_L tan^/Zcy ,-65 > 1 if the pitch angle ^ is not close to 0, we can simply use the 
asymptotic form of Bessel function for large argument. Then we can derive an analytical 
result for Alfvenic turbulence (see Appendix B) , 



D 
D 



pp 



,2.5 



cos ^ 



5.5 



^1.5^2.5 sin ^ 



T[6.5,k 



maxkresL^l 



(10) 



where r[a, z] is the incomplete gamma function. The presence of this gamma function in 
our solution makes our results orders of magnitude larger than those^ in Chandran (2000) 
for the most of energies considered. However, the scattering frequency. 



iy^2D,,/{l-fx'), (11) 

are much smaller than the estimates for isotropic and slab model (see Yan & Lazarian 2002). 
The Alfven modes are damped in partially ionized medium at scales larger than the resonant 
scales of moderate energy CRs. As the anisotropy of the Alfven modes is increasing with 
the decrease of scales, the interaction with Alfven modes becomes more efficient for higher 
energy cosmic rays. When the Larmor radius of the particle becomes comparable to the 
injection scale, which is likely to be true in the shock region as well as for very high energy 
cosmic rays in diffuse ISM, Alfven modes get important. 

It's worthwhile to mention the imbalanced cascade of Alfven modes (CLV02). Our basic 
assumption above was that Alfven modes were balanced, meaning that the energy of modes 
propagating one way was equal to that in opposite direction. In reality, many turbulence 
sources arc localized so that the modes leaving the sources are more energetic than those 
coming toward the sources. The energy transfer in the imbalanced cascade occurs at a slower 
rate, and the Alfven modes behave more like waves. The scattering by these Alfven modes 
is hkely to be more efficient. However, as the degree of anisotropy of imbalanced cascade is 
currently uncertain, and the process will be discussed elsewhere. 



^We compared our result with the resonant term as the nonresonant term is spurious as noted by Chandran 
(2000). 
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4.2. Scattering by fast modes 

The contribution from slow modes is no more than that by Alfven modes since the 
slow modes have the similar anisotropics and scalings. More promising are fast modes, 
which arc isotropic (CL02). With fast modes there can be two types of resonant interaction: 
gyroresonance and transit-time damping (TTD) (Schlickeiser & Miller 1998). 

Fast modes potentially can scatter CRs by transit-time damping. From the resonant 
condition uo — k\\V\\ ~ 0, we see that the contribution is mostly from nearly perpendicular 
propagating modes (cos^ ~ 0). According to Eq.(4),we see that the corresponding cor- 
relation tensor for the magnetic fluctuations are very small, so the contribution from 
TTD to pitch angle scattering is less important unless the pitch angle is close to 90°. The 
advantage of TTD is that it doesn't have a distinct resonant scale associated with it. The 
resonant condition only requires k\\/k = Vph/V\\. As the result, TTD is an alternative to 
scatter low energy CRs whose Larmor radii are below the damping scale of the fast modes. 
Moreover, TTD can be substantial to the momentum diffusion or acceleration (also known 
as the second order Fermi acceleration). It can be crucial in some circumstances, e.g., for 
7 ray burst (Lazarian et al. 2003), and super- Alfvenic grain acceleration (Yan & Lazarian 
2003). 

Here we apply our analysis to the various phases in ISM. 

Recent observations (Beck 2001) suggest that the galactic halos are magnetic-dominant, 
corresponding to low f3 medium. For low /3 medium, we use the tensors given by Eq.(4). For 
the halo, the collisionless damping is dominant (see Table 1). As we see from Eq.(A2), the 
damping increases with 6 unless 6 is close to 7r/2, which results in kj_^c < res- This means 
that the fast modes propagating at smaller angle with the magnetic field are less damped. 
As a result the argument for the Bessel function in Eq.(Bl) is k± tan ^/k^^^res < 1 unless ^ is 
close to 90°. So we can take advantage of the anisotropy of the damped fast modes and use 
the asymptotic form of Bessel function for small argument J„(x) ~ (x/2)"/n! to obtain the 
corresponding analytical result for this case (see Appendix B): 



D 



0.5/ 



D. 



pp 



2m 



'I - 



.res 

2 



+ l)-^)/7 



(12) 



Prom this equation, we see the key factor is {k±^c/ k\\,res) — tan^c- As the coUisionless 
damping is a function of 9, the cutoff scale depends also on 9. Put the Eq.(A2) into the 
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relation TfeFi ~ 1, we see cos}-^ 9c/ s\v? 9c oc k^^^l^g oc 7"°-^. Therefore scattering frequency 
in the halo increases slowly with energy as shown in Fig. la. For TTD, the contribution is 
mostly from the nearly perpendicular propagating modes for which the coUisionless damping 
is small (see Eq.A2). Thus the asymptotics of the Bessel function is not applicable. We 
obtain a numerical solution from Eq.(B3). The contribution of TTD to scattering is much 
smaller than that of gyroresonance. 

Hot ionized medium (HIM) is in high (3 regime and we adopt the tensors given by 
Eq.(8). The mean free path of the HIM is of the order of parsec and thus the fast modes 
in this medium are subjected to coUisionless damping (sec tablel). We use Eq.(A4) to get 
the corresponding damping scale k^^ = 1.3 x 10^^ sin^^'^ ^/ cos^^cm. Therefore, fast modes 
become also anisotropic in HIM resulting in A:^ < fcy ^es- For CRs with energy Ej. < 25GeV, 
k\ires > k±^c = \Jk1 — /cy , SO wc Can use the asymptotics of Bessel function J„(a;) ~ (a;/2)"/n! 
to get an analytical result for gyroresonance in high [3 medium: 



27r(QT;/x)°-5(l-At^) 

/I. \7/2 /. \ 11/2 

/t. \3/2 /i. \7/2 ' \^^) 

)mVi/3/3 - (1 - (^) )m^Vl/7(5 _ 

from which we see the scattering frequency decreases with k\\^res a-nd so increases with Ej^. 
This result agrees well with the numerical solution obtained by Eq.(9) (see Fig. la). At high 
energy end, the resonant scales are much larger than the damping scale k~^ so the effect 
of damping to the scattering curve can be ignored. As a result, the shape of the scattering 
curve is similar to the case where there is no damping (see Yan & Lazarian 2002). The 
scattering frequency in this case according to Eq.(9), increases with Larmor frequency f2o/7 
and the resonant scale k~}g ~ , which increase with proton energy Ek- For Ek > IGeV, 
7 goes nearly linearly with E^ so that the scattering frequency decreases due to the flo/j 
dependence which changes more rapidly. 

According to the results in §2, the warm ionized medium (WIM) is also in low f3 but 
collisional regime. The mean free path Imfp = 6 x lO^^cm (tablel), which corresponds to the 
resonant scale of CRs with energy ~ lOGeV. Thus For CRs harder than lOGeV the resonant 
modes with k\\ — ^/v\\ are subjected to collisional damping. Since the viscous damping 
increases with 9, we can apply Eq.(12) to the gyroresonance in WIM as well. As the viscous 
damping in the low P medium scales as /c^. while the cascading rate t^^ oc k^'^, we can get 
cos-*^-^ 9c/ sin^ 9c oc kl'^^^ oc 7"^-^ from the truncation condition TfcFi ~ 1. Thus according to 
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Eq.(12) the scattering frequency curve in WIM gets steeper (see Fig. la). For lower energy 
CRs, the only available modes are those with k in a small cone as the residual of the coUisional 
damping at large scales. As a conservative estimate, wc adopt the smallest perpendicular 
viscous cutoff scale kj_^c above the mean free path to do the calculation. Certainly, the 
resonant modes corresponding to this energy range of CRs are also affected by collisionlcss 
damping. However, the pitch angle of these modes are so small that the coUisionless process 
is marginal (see Eq.(A2)). 

In partially ionized medium, the fast modes arc severely damped by the ion-neutral 
damping. The cascade is cut off before the resonant scales for most of the energy range we 
considered. Indeed, for the parameters chosen, we find that gyrorcsonancc only contribute 
to the transport of CRs of energy > O.STeV in DC and > ITeV for WNM and CNM. 
Therefore only TTD contributes to the scattering of moderate energy CRs in these media. 
Since k±v±/ft < kctsin ^ / k\\^res < 1, we can use the asymptotics of Bessel function J„(x) ~ 
(x/2)"/n! to get an analytical result for TTD in these media (see Appendix B) 
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(14) 



from which we see D^^ is approximately ccl/v. Ks v increases with E^; when E^. < rUpC^ ~ 
IGeV and ~ c when E^ > IGeV, the scattering frequency decreases with energy when 
Ek < IGeV and keeps nearly a constant at higher energies. As expected, the scattering is 
much less efficient in partially ionized media because of ion-neutral damping. 

From the above results, we see that the fast modes dominate CRs' scattering in spite 
of the damping. From Eq.(14), we see the ratio of scattering and momentum diffusion 
rates p^D^^/ Dpp — 1///^, which means that the scattering and acceleration by TTD are 
comparable. However, for gyroresonance, we know that p^D^^/Dpp ~ v'^/{min{f3,l}V2) 
(see Eq.(13)and Eq.(12)). Consequently for momentum diffusion (or acceleration) TTD is 
dominant and the rate differs from the scattering rate by a factor ~ /x^. 

A special case is that the cosmic rays propagate nearly perpendicular to the magnetic 
field, so called the 90° scattering problem. The dynamics of the turbulence was suggested to 
be taken into account in this case (Schlickciscr & Achatz 1993). There have been attempts 
of using the concept of the dynamical turbulence in dealing with the CR transport before 
(Bieber et al. 1994). It has been shown in Schlickeiser & Achatz (1993) that if the decay 
time akVA is longer, the resonance function can be approximated by 5 function. In other 
words, et ^ is the "magnetostatic hmit". From (Eq.(3)), we can get a — (/cL)"^/^. 
For moderate energy CRs, the corresponding a « 1 so that S function is still a good 
approximation even for the 90" case. The corresponding resonant scale is ~ Vph/Q. Thus 
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unlike with Alfven modes, we can implement unified 5 function and then integrate over the 
pitch angle /x to estimate the spatial diffusion. Here we focus on galactic halo since the 
scattering rate is the highest there^. In the precedent calculations, we ignored Vph because 
it is negligible comparing with v/x. Now should be replaced with vfi ± Vph as seen from 
the gyroresonance condition u — k\\V\\ = nQ. The complication comes from the cutoff due to 
damping. As addressed earlier, the fast modes become slab-type at small scales. Thus we 
can get 6c from Eq.(3) and Eq.(A4) k±^c/k\ires — Oc — 2{miV\\/ nrriefiVLlY/^ ^ 1. Then we 
can obtain an analytical solution from Eq.(12) for the parallel spatial diffusion coefficient: 



1 -^MM 

3/-?; / mg/?^ 



t(;^) [-l+2. + 2(-l + ain^], (15) 

where e = Va/v. This result is compared with numerical integral of Eq.(12) in Fig. 2b. These 
results are not applicable to high energy CRs for which the slab-model of fast modes fails. 

It should be pointed out that all the results above are very much dependent on the 
plasma /3, which is somewhat uncertain. But the basic conclusion is definite: fast modes 
dominate the cosmic ray scattering except in very high /5 medium where fast modes are 
severely damped. This demands a substantial revision of cosmic ray acceleration/propagation 
theories, and many related problems should be revisited. 



5. Cosmic ray self confinement by streaming instability 

When cosmic rays stream at a velocity much larger than Alfven velocity, they can 
excite resonant MHD modes which in turn scatter themselves. This is so called streaming 
instability. It was claimed that the instability could provide confinement for cosmic rays with 
energy less than lOOGeV (Cesarsky 1980). However, this was calculated in an ideal regime, 
namely, there was no background MHD turbulence. Only ambipolar coUisional damping was 
involved. In other words, the self-excited modes would not be damped in fully ionized gas, 
which is not true for turbulent medium. Here we shall reconsider in turbulent ISM how this 
instability would affect cosmic ray transport, whether it is still an effective mechanism or 
not. 



^It looks from Fig. la that the scattering is most efficient in WIM for CR>10GeV. However, it's just for 
a particular ^i. For small /x, which corresponds to small scales and is the most important for the calculation 
of mean free path, the scattering in halo is still dominant. 
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Fig. 1. — Scattering by fast modes. The scattering frequency u given by Eq.(ll) vs. (a) the 
kinetic energy of cosmic rays, (b) kinetic energy E]^ of cosmic rays in partiaUy ionized 
medium. The sohd hues represent our analytical results. In (a), the diamond line is our 
numerical result for galactic halo, the 'o' line refers to HIM and the 'x' line represents WIM. 
In (b), the diamond line is our numerical result for WNM, the 'o' line refers to CNM and 
the 'x' line represents MC. 




Fig. 2. — (a)The scattering frequency v vs. cosine of pitch angle ji of IGeV CR in fully 
ionized medium. The solid lines represent our analytical results, the diamond line is our 
numerical result for galactic halo, the 'o' line refers to HIM and the 'x' line represents WIM. 
(b). The mean free path of CRs in galactic halo resulting from numerical integration of 
Eq.(12) 
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The growth rate of the modes of wave number k is given by (Longair 1997) 



r{k) = u 



N{> E) 




(16a) 



'0 



n, 



'p 



where A'"(> E) is the number density of cosmic rays with energy > E which resonate with 
the wave, i.e., v{E) = ^l/k\\, Hp is the number density of charged particles in the medium. 
If taking the energy spectrum of cosmic ray particles n{E) ~ 1.6 x 10^(£'/MeV)~^''^ as 
measured at the Earth (Longair 1997), then A^(> E) ~ 0.25 x 10-^°{E/GeV)-^-'^cm-^sr-^. 

The damping rate due to the turbulence cascade depends on the properties of turbulence. 
To estimate the rate at which background turbulence can suppress th streaming instability we 
consider both interaction with fast and Alfvenic modes. The interaction with the fast modes 
happens at the rate ~ {k/ L)~^/^Vph/V'^ (see Eq.(3)). The interaction with Alfven modes 
is less straightforward to quantify because of their anisotropy. Consider the interaction of a 
perturbation created by streaming particles and the background Alfvenic turbulence. If we 
approximate this perturbation by a wave with k — k\i, it gets evident that the interaction 
will be very different from what takes place between the oppositely moving wave packets in 
Alfvenic tTirbulcncc. If we consider local in k-spacc interaction the plane wave created by the 
streaming instability will undergo random shearing. The mean variation of the phase A(/9 will 
be approximately 1/\/N, where ~ {k±/k\\y. The decorrelation time tdc can be estimated 
using random walk arguments (A(^)^ • (t^c'^) ~ 1, which provides tdc ^ i^~^{k±/k\\y. 

Thus the dominant damping mechanism here is the cascade of fast modes. By equating 
the growth rate Eq.((16a)) and the damping rate (Eq.(3)), we can find that the streaming in- 
stability is only applicable for particles with energy < 1.5xlO~^[n~^{v/VAY'^{L^lo/VA)^'^]^^^''^GeV. 
In HIM, the corresponding energy is ~ SOOGeV. As shown in the Eq.(16a) , the growth rate 
depends on the CRs' density. In those regions where high energy particles are produced, 
e.g., shock fronts in ISM, 7 ray burst, SN, the streaming instability is more important. 



We have attempted to apply a realistic model recently obtained for MHD turbulence 
to the problem of cosmic ray scattering taking into account the turbulence cutoff arising 
from both coUisional and coUisionless damping. We considered both gyroresonance and 
transit-time damping (TTD), and used QLT to obtain our estimates. 

In earlier papers it was frequently assumed that Alfvenic turbulence follows Kolmogorov 
spectrum while the problem of CR scattering was considered. This is incorrect because 



6. Discussion 
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MHD turbulence is different from hydrodynamic turbulence if magnetic field is dynamically 
important. Realistic MHD turbulence can be decomposed into, Alfven, slow and fast modes. 
They have different statistics (see CL02). The Alfven and slow modes follow GS95 scalings 
and show scale-dependent anisotropy. The anisotropy of Alfven and slow modes reduce the 
scattering efficiency to such a degree that their impact on cosmic ray propagation is marginal. 

Fast modes, however, are promising as a means of scattering due to their isotropy. Even 
in spite of the damping fast modes dominate scattering if turbulent energy is injected at large 
scales. We provided calculations for various phases of ISM, including galactic halo, HIM, 
WIM and partially ionized media. As the fast modes are subjected to different damping 
processes, the CR scattering varies . In galactic halo and HIM, the fast modes are subjected 
to coUisionless damping. The mean free path in WIM is much smaller, and ion viscosity 
dominates the damping. All these dampings decrease with 6 and marginal for parallel modes. 
As the result, the fast modes at small scales become anisotropic with k\\ » k±. In this sense, 
fast modes are in slab geometry but with less energies. These modes can efficiently scatter 
CRs with the resonant scale k\\^res = ^/'^\\ larger than the cutoff through gyroresonance. 
TTD is usually weaker for pitch angle scattering than gyroresonance, but dominant for 
momentum diffusion. In all the cases, the scattering efficiency decreases with plasma j3 
because the damping increases with j3 and the perturbations of turbulence decreases with (3. 

How good is our model of ISM turbulence? We used spectra for Alfven and fast modes 
which follow from theoretical constructions and which were tested numerically (CL02). 
Within the ISM similar spectra develop dynamically when random driving is applied at 
large scales (see Mac Low 2003) and therefore we expect that the distribution of energy over 
scales for different modes to coincide with the assumed ones. The partition of energy over 
modes depend, however, on driving. For random driving with 5V ~ Va we expect approxi- 
mately equal energy in fast and Alfven modes. As the details of ISM driving are not clear 
and the QLT provides results correct up to a factor of the order unity, our predictions for fast 
mode scattering contain an uncertain factor of the order unity. We believe, however, that our 
results for scattering by Alfven modes that given by Eq.(lO) are subject to more uncertain- 
ties related to the turbulence driving. Indeed, the scattering rates that follow from Eq.(lO), 
although orders of magnitude larger than those in Chandran (2000), are still extremely small 
compared to models with isotropic spectra. Because it is so small, higher order effects may 
become important. Does the steady state anisotropic spectrum has always enough time to 
develop? The dynamic establishing of anisotropic spectra of Alfvenic turbulence requires a 
few eddy turnover times. For imbalanced driving this can take somewhat longer time (see 
discussion in §5). Therefore one can imagine situations when the spectrum deviates of its 
steady state form given by Eq. (3). These deviations are most likely to be small. There- 
fore we do not expect that in typical ISM conditions the contribution of Alfvenic modes is 
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comparable to that of the fast modes. Identification of the particular situations when the 
opposite may be true requires a more systematic study of particular driving mechanisms and 
is the subject of future research. We also note, that if energy is injected at the small scales, 
the anisotropy of the Alfvenic modes is not large and the scattering can be appreciable. The 
deviations of the spectrum from the steady state one are not important. 

Although scattering via fast modes is isotropic at large scales and becomes slab-type at 
small scales at which damping gets important, our model is radically different from earlier 
discussed slab and isotropic+slab models (Bieber et al. 1994, Schlickeiser & Miller 1998). We 
used theoretically motivated models of turbulence that have been tested numerically. There- 
fore the transitions from one regime to another as well as the values of the intensity of 
perturbations have good justification within our approach. 

We adopted a particular set of parameters of ISM while doing the calculation. The 
basic conclusion, namely, fast modes dominate the cosmic ray scattering, will remain the 
same if the parameters are altered. However, since the (5 value in ISM remains uncertain 
from observation, the damping scale kc is somewhat uncertain. And this determines the 
energy limit down to which gyroresonance is dominant. 

7. Summciry 

In the paper above we characterized interaction of cosmic rays with balanced interstellar 
turbulence driven at a large scale. Our results can be summarized as follows: 

1. Fast modes provide the dominant contribution to cosmic ray scattering. 

2. Scattering of cosmic rays varies from one interstellar phase to another. It sensitively 
depends on the damping of fast modes. 

3. Due to anisotropic plasma damping fast modes develop at small scales anisotropy 
which makes gyroresonant scattering within slab approximation applicable. At scales where 
damping is negligible the isotropic gyroresonant scattering approximation is applicable. 

4. In partially ionized gas fast modes scatter cosmic rays of energies < ITeV by TTD 
mechanism. 

5. Streaming instability is partially suppressed due to the interaction of the emerging 
magnetic perturbations with the surrounding turbulence. 

We acknowledge valuable discussions with Jungyeon Cho, Steve Cowley, Don Cox, 
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Below we summarize the damping processes that we consider in the paper. 

Ion-neutral damping 

In partially ionized medium, viscosity of neutrals provides damping (see LY02). If the 
mean free path for a neutral atom, /„, in a partially ionized gas with density Utot = + Ui 
is much less than the size of the eddies under consideration, i.e. Ink <C 1, the damping time 



where z/„ is effective viscosity produced by neutrals^. The mean free path of a neutral atom 
In is influenced both by collisions with neutrals and with ions. The rate at which neutrals 
collide with ions is proportional to the density of ions, while the rate at which neutrals 
collide with other neutrals is proportional to the density of neutrals. The drag coefficient for 
neutral- neutral collisions is ~ 1.7 x 10~^°r(K)°-^ cm^ s~^ (Spitzer 1978), while for neutral- 
ion coUisions it is ~ {vrOin) ~ 1.9 x 10~^ cm^ s~^ (Draine, Roberge & Dalgarno 1983). 
Thus collisions with other neutrals dominate for Ui/rin less than ~ 0.09T°-^. Turbulent 
motions cascade down till the cascading time is of the order of tdamp- The maximal damping 
corresponds to k^-^ ~ Z„. If the neutrals constitute less than approximately 5%, the cascade 
goes below and is damped at smaller scale (see below). 

Collisionless damping 

The nature of coUisionless damping is closely related to the radiation of charged particles 
in magnetic field. Since the charged particles can emit plasma modes through acceleration 

(cyclotron radiation) and Cherenkov effect, they also absorb the radiation under the same 
condition (Ginzburg 1961). The thermal particles can be accelerated cither by the paral- 
lel electric field which can also be called Landau damping or the magnetic mirror (TTD) 
associated with the co moving modes (or under the Cherenkov condition A;||f|| = u). The 
gyroresonance with thermal ions also causes the damping of those modes with frequency 



A. A. Damping of MHD turbulence 




(Al) 



^The viscosity due to ion-ion collisions is typically small as ion motions are constrained by the magnetic 
field. 
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close to the ion-cyclotron frequency (Leamon et al. 1998), though it is irrelevant to the 
low frequency modes since we deal with GeV CRs here. The coUisionless damping depends 
on the plasma (3 = Pgas/Pmag and the propagation direction of the modes. In general, the 
damping increases with the plasma j3. And the damping is much more severe for fast modes 
than for Alfven modes. For instance, the Alfvcn modes arc weakly damped even in a high 
(3 medium, where fast modes are strongly damped. The damping rate 7^ = r^^ of the fast 
modes of frequency a; for /? <S 1 and ^ ~ 1 (Ginzburg 1961) is 



+ 5 exp( 



4 cos 6* 
1 



rriH 



exp( 



tuhP cos^ 9 



) 



)], 



(A2) 
(A3) 



where mg is the electron mass. The exact expression for the damping of fast modes at small 
9 was obtained in Stepanov® (1958) 
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When /3 ^ 1, we obtain the damping rate as a function of 9 from Foote & Kulsrud 1979, 



^ sin^^ r 2uj^/Q, for k < n^/ (5Va 
^ cos3 9^\ 20^7/3 for k > (5Va ^ ^ 

where is the ion gyrofrequency. 

Ion viscosity 

In a strong magnetic field {VtiTi ^ 1) the transport of transverse momentum is prohib- 
ited by the magnetic field. Thus transverse viscosity r^j, is much smaller than longitudinal 
viscosity rjo = 0.96n(T/eV)Ti, rjx_ ~ rjo / {VtiTiY' . The heat generated by this damping is 
Qvis = Vo{9vx/dx + dvy/dy — 2dvz/dzy /3, where v^, Vy, are the velocity components 
of the wave perturbation (Braginskii 1965). From the expression, we see that the viscous 
damping is not important unless there is compression. Therefore Alfven modes is marginally 
affected by the ion viscosity. 



We corrected a typo in the corresponding expression. 
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While the damping due to compression along the magnetic fields can be easily under- 
stood, it is counterintuitive that the compression perpendicular to the magnetic also results 
in damping through longitudinal viscosity. Following Braginskii (1965), the damping of per- 
pendicular motion may be illustrated in the following way. To understand the physics of 
such a damping, consider fast modes in low f3 medium. In this case, the motions are pri- 
marily perpendicular to the magnetic field so that dvx/dx = h/n ~ B/B. The transverse 
energy of the ions increases because of the conservation of adiabatic invariant v\/B. If the 
rate of compression is faster than that of collisions, the ion distribution in the momentum 
space is bound to be distorted from the Maxwellian isotropic sphere to an oblate spheroid 
with the long axis perpendicular to the magnetic field. As a result, the transverse pres- 
sure gets greater than the longitudinal pressure by a factor ~ n/n, resulting in a stress 
~ pTih/n ~ rjoOvx/dx. The restoration of the equilibrium increases the entropy and causes 
the dissipation of energy with a damping rate Fjo^ = kj_r]Q/6pi (Braginskii 1965). 

In high f3 medium, the velocity perturbations are radial as pointed out in §2. Thiis 
Qvis — VoikxVx + kyVy — 2kzVzY/'i = i]ok'^v'^{sm'^ 6 — 2cos^^)^/3. Dividing this by the 
total energy associated with the fast modes Ek — Piv"^, we can obtain the damping rate 

All in all. 



P ^ f klr]o/Qpi for/5<l 

\/t2r?o(l-3cos2^)V6p, for/3>l ^ ^ 

Resistive damping 

In the paper we ignore the resistive damping because of the following reason. The 
resistivity is much smaller than the longitudinal ion viscosity for fast modes. For parallel 
propagating fast modes which are not subjected to the viscous damping from the longi- 
tudinal viscosity, the resistive damping scale can be obtained by equating the damping 
rate c^k'^/87ia± with the cascading rate of the fast modes, where a± ~ 10^^(T/eV)'^/^ (see 
Kulsrud & Pearce 1969) is the conductivity perpendicular to the magnetic field: k~Jg = 
1.8 X 10'^(L/pcy^'^(T/10^K)~^{VA/kms~^)~'^^^cm. This scale turns out to be much smaller 
than the mean free path, where collisionless damping takes over. For Alfven modes, the 
resistive damping rate is Tres = c^k'^/8na\\ (the term associated with /cyis negligible be- 
cause of anisotropy), where a\\ — 1.96(7xis the conductivity parallel to the magnetic field. 
By equating it with the cascading rate of Alfven modes, we can get the resistive damp- 
ing scale k~Jg ~ L(c^/87r(T||VA-^)^'^^- Comparing with the proton Larmor radius, we find 
k-\/rL = 3.6 X 10-3(L/pc)V4(r/10^X)-=^/2^V2(y^/Cji/4 j^'g clear that in many astro- 
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physical plasmas the resistive scale is less than proton Larmor radius, where Alfven modes 
can not proceed further because of anomalous resistivity. All in all, we see that the damping 
arising from resistivity is marginal. 



B. B. Fokker-Planck coefficients 

In quasi-hnear theory (QLT), the effect of MHD modes is studied by calculating the 
first order corrections to the particle orbit in the uniform magnetic field, and the ensemble- 
averaging over the statistical properties of the MHD modes (Jokipii 1966, Schlickeiser & 
Miller 1998). Obtained by applying the QLT to the coUisionless Boltzmann-Vlasov equation, 
the Fokker-Planck equation is generally used to describe the evolvement of the gyrophase- 
avcragcd particle distribution. The Fokker-Planck coefficients -D^^, -D^p, Dpp arc the funda- 
mental physical parameter for measuring the stochastic interactions, which are determined 
by the electromagnetic fluctuations: 

< 5„(k)5;(k') >= 5(k - k')Pa^(k) < 5„(k)^^(k') >= 5(k - k')T«/3(k) 

< E«(k)5;(k') >= 5(k - kOga^(k) < E,(k)^;(kO >= (5(k - k')i?a/3(k), 

Adopting the approach in Schlickeiser & Achatz (1993), we can get the Fokker-Planck 
coefficients. 
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where |kmin| = ^mm = L~^, Ik^axl = kmax corresponds to the dissipation scale, where 
C = 1 for Alfven modes, dk^ = 27[k±dk±dk\\, and ( = k\\/k for fast modes, m = jmn is 
the relativistic mass of the proton, v± is the particle's velocity component perpendicular to 
Bo, (f) = a,rcta,n(ky / kx) , C,TZ = (x ± iy)/^/2 represent left and right hand polarization. For 
low frequency MHD modes, we have from Ohm's Law E(k) = — (l/c)v(k) x Bq. So we 
can express the electromagnetic fluctuations Tjj, Rij in terms of correlation tensors Kij. The 
particular dispersion relations for cu are not important unless scattering happens at an angle 
close to 90°. For uj we use dispersion relations u; — k\\VA for Alfven modes and u — kVph for 
fast modes. Those dispersion relations were used to decompose and study the evolutions of 
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MHD modes in CL02. 

For gyroresonance, the dominant interaction is the resonance at n = ±1. By integrating 
over r, we obtain 
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For TTD, 
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Define the integrands 



r(k) = 2Jl{^){{Knn{\^)-rKcc{\^)-r2U^)U^){e^^^K^^^ (B4 

. We shall show below how they can be simplified in various cases to enable an analytical 
evaluation of the integral. The spherical components of the correlation tensors are obtained 
in the following. 
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For Alfven modes, their tensors are proportional to 



^ . sin^ (f) — cos (f) sin ( 
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and 



^'''Inc + e-'^^cn = ^^^^^ x e^^^ + ii^L±!^ x e"'^^ = - cos 20 + + sin 20 = - 

For Alfven modes, k±v±/fl 3> 1 because of the anisotropy. So we can use the asymptotics 
of the Bessel function J„(a;) = ^y2/-KX cos(x — n7r/2 — 7r/4) when x oo. Thus from 
Eq.(6,B5,B5), we simplify the integrand in Eq.(9), G(k) oc cxp(— L-'^/^|A:|| |/A:^^"^), the 

integral of which will be oc {\k\\\L'^^^)~^'^T[Q.5, kmax 

Tj^g^ from Eq.(9) we can get 
the analytical result for the scattering by Alfven modes as given in Eq.(lO). 

For fast modes, their tensors have such a component 



H,,^Ak-'-'( '""f^ ^^^'^f . (B5) 
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Thus we have 
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and 



(^''^Hnc + e-'''^Hcn = (^x. - Hyy) cos 20 + {H,y + Hy,) sin 20 = 1. (B7) 
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In general, it's more difficult to solve the integral in Eq.(Bl) for the fast modes because 
they are isotropic. However, if taking into account damping, Bessel function can be evaluated 
using the zeroth order approximation. Thus from Eq.(9,B6,B7), we see for gyroresonance 
the integrand can be simplified as G(k) oc (J2 ~ JoYk~'^^'^ oc A;"^/^. For TTD, T(k) oc 
Jf{k±vj_/Q)k^'^^'^ oc k^^^'^v'j_sm'^ 6 oc k~^^'^vj_{l — V^f^/v'^^). Then we can solve the integral in 
Eq.(B3) analytically and get the scattering rate for fast modes. 
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